Generalized Fermat principle for classical and quantum systems 
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The analogy between dynamics and optics had a great influence on the development of the foun- 
dations of classical and quantum mechanics. We take this analogy one step further and propose a 
new variational principle in classical phase space and quantum Hilbert space that mimics Fermat 
principle in optics. The state of a many-body classical (quantum) system follows a trajectory in the 
phase (Hilbert) space that extremizes the time functional. The classical Hamilton's equations of 
motion and the quantum Schrodinger's equation are shown to lead necessarily to a stationary time 
functional. A plethora of other unconventional variational principles related to the geometry of the 
dynamical evolution, and inspired from Fermat principle are discussed in this paper. 



It has been a hallmark in the history of physics in the 
past few centuries that ingenious and elegant reconsid- 
erations of the relations between geometric concepts and 
physical phenomena lead to paradigm shifts in physics 
and our perception of the universe. The most notable ex- 
amples are the theories of special and general relativity. 
Novel reconsiderations of geometric quantities in physics 
often lead to new variational principles which assign the 
natural evolution of physical systems with an extremum 
of some functional or a geodesic curve in some hyper- 
space. 

The quest to find new variational principles in physics 
has not stopped till today pQ. The oldest of these vari- 
ational principles is Fermat principle of least time which 
became a fundamental principle in geometric optics. The 
principle was introduced by Fermat, which he also called 
the principle of natural economy [2] , and states that light 
rays travel in a general medium along the path that min- 
imizes the time taken to travel between the initial and 
final destinations. The concept of natural economy in- 
spired Maupertuis to introduce the principle of least ac- 
tion in analytical mechanics which later evolved through 
the work of Euler, Lagrange, Hamilton, and Jacobi to 
become a fundamental and universal concept in classical 
mechanics. 

The mathematical formulation of Fermat principle 
states that the time functional T defined as 

r-/"-£ (i) 

is minimized along the trajectory of light rays [3] . In other 
words, if T is computed along all possible trajectories be- 
tween an initial and final positions Sf where v(s) is the 
speed of light and ds is the distance element along the 
trajectory, T will always be minimum along the actual 
path travelled by the light rays. The modern version of 
Fermat principle is written in terms the index of refrac- 
tion n(s) = where c is the speed of light in free 

space and states that the optical path length n(s)ds 
is a minimum. In that sense, Fermat principle is the 
optical analog of Jacobi least action principle [4], which 



states that for a conservative classical system at energy 
E, with potential function between its constituent parti- 
cles V, the action functional 



1=1 - V(s)ds (2) 
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computed in the configuration space along the physical 
path traversed by the system between the initial and fi- 
nal representative points is an extremum. The remark- 
able property of this action that distinguishes it from 
other variational principles in analytical mechanics, i.e., 
Hamilton's and Lagrange's variational principles, is that 
it represents a purely geometric quantity that is com- 
puted along different trajectories in the configuration 
space between fixed points without referring to any evolu- 
tion in time. Therefore, equation|2]can be used to define a 
new Riemannian space whose metric ds' = \/ E — V(s)ds 
where the natural evolution of the representative point 
of the system will be along geodesic curves. 

In this work we elevate the status of the time func- 
tional to a universal action in both classical and quan- 
tum generic many-body systems that yields an extremum 
along the physical path in a proper space S that fully 
characterizes the state of the system and its future evo- 
lution, once the Hamiltonian is known. A proper space 
for quantum systems is the projective Hilbert space P, 
where wavefunctions are defined upto an overall phase 
factor, and for classical systems a proper space is the 
phase space consisting of all the generalized coordinates 
and momenta {qi,Pi}. The principle informally states 
that among all possible trajectories between Si and Sf 
only those which make T invariant under an infinitesi- 
mal variation in the path, i.e., 

ST = (3) 

are possible candidates for the dynamical evolution. 

The expression for the speed of the representative point 
of the system v(s) computed over S is obtained from the 
equations of motion of the system that will drive its evo- 
lution along the physical path (i.e., Schrodinger's equa- 
tion in quantum systems and Hamilton's equations of 
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motion in classical systems). We emphasize here that we 
are not on the mission of deriving these equations of mo- 
tion from the time action, because we have to use them 
to find v(s) in the first place, but rather to assert that 
they necessarily lead to a stationary time action. Despite 
being a purely geometric quantity, the time functional 
by definition can straightforwardly extract the timing in- 
formation of the evolution once the path is determined. 
However, Fermat principle in a generic many-body sys- 
tem is in general defined for a time-independent system, 
i.e., the Hamiltonian does not contain time explicitly. 
Nevertheless, we use insights gained from this principle to 
define a new variational principle for quantum mechanics 
which is valid for a generic time-dependent Hamiltonian. 
Our study opens the door for introducing new variational 
principles in other classes of physical systems not consid- 
ered in this paper (e.g., continuum mechanics and field 
theory) and hopefully shed more light on the role of time 
in physics. 



— Fermat Principle in Hilbert Space Unlike the the- 
ory of general relativity, where geometric concepts are 
rooted in the bones of the theory, quantum mechanics 
did not enjoy such a solid connection with geometry in 
its very beginning. Later development of the geometric 
phase in quantum mechanics has encouraged many re- 
searchers to pay more attention to this issue. Anandan 
and Aharonov have initiated the investigations of the na- 
ture of the geometry of quantum evolution in the projec- 
tive Hilbert space P through a series of papers in the late 
80s [5HZ]- They have shown [6] that the speed of quan- 
tum evolution in P is related to the energy uncertainty 

AE= ((H 2 )-(H) 2 y /2 via 



ds = 2AE dt/h 



(4) 



where ds is the infinitesimal distance in P given by the 
Fubini-Study (FS) metric ds 2 = (5^\5^) - | (<fy#) | 2 . 
The trajectory traversed by a ray in P under unitary 
evolution is not a geodesic, i.e., S j ds ^ 0. This can be 
immediately envisaged by considering a system composed 
of a single quantum spin 1/2. In this case, P is simply 
the Bloch sphere and the precession motion of the spin 
on Bloch sphere is not a geodesic. Several attempts [9] 
tried to conceive new formulations where the quantum 
evolution is a geodesic flow. 

The simple answer to this problem suggested by the 
present paper is to consider T = J as a geodesic 
quantity, which is stationary along the physical path 
through P, i.e., Fermat Principle in Hilbert Space. This 
is to be distinguished from the quantum brachistochrone 
problem [TO], where the Hamiltonian that leads to opti- 
mal time evolution between an initial and final state is 
sought. The above proposition states that the unitary 
evolution generated by any Hamiltonian is optimal, with 



respect to all other possible trajectories connecting the 
initial and final states (Fig. l.a). 

If we parametrize the evolution along any path con- 
necting and \ipf) by some arbitrary parameter r, we 
can write equation [T] as 

T- t'\r , (5) 
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where \ip) = and P is the projection operator l^)^. 
Taking the variational derivative of T with respect to 
(5ip\ subject to the constraints of normalization and fixed 
initial and final states, we arrive at the Euler-Lagrange 
(EL) equation of the integrand L, 
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After taking the inner product with |^) and making use 
of the constraint condition (^1^) = 1 we arrive finally at 



(m 2 = 



<V#)-|<V#>| 2 



mf (7) 



We can see straightforwardly that the Schrodinger's 
equation \ip) = ±iH\ip) is the only linear evolution for 
\ift(t)) that will satisfy this condition. The fact that both 
the forward and backward time evolutions exist simul- 
taneously on the energy shell for the same variational 
problem, i.e., between the same initial and final states, 
may be related to the quantum ergodic theorem. Systems 
possessing the quantum ergodicity property as defined in 
[TT] are those systems "where all states, within a given 
energy range can be reached from all other states within 
the range even if only via other intermediate states" . In 
our case, if quantum ergodicity holds, it will be possi- 
ble to travel from any |^) to any \ipf) via two different 
routes. 

Substituting Schrodinger's equation directly into EL 
equation would have been the easier way to show the va- 
lidity of the principle. However, we highlight equation [7] 
to point out the intriguing possibility that some nonlin- 
ear equations may satisfy this equation and consequently 
fulfills Fermat principle. This is indeed the case, and we 
found that the nonlinear equation 



(8) 



satisfies equation [7| 

As emphasized above, the timeless nature of this vari- 
ational principle does not make it adequate to deal with 
time-dependent Hamiltonians. To achieve this aim, we 
need to add a 'time-like' component to the Fubini-Study 
metric and modify the velocity term v(s) accordingly. 
There exists a variety of ways to implement this pre- 
scription. We adopt the following choice of the full met- 
ric, ds' 2 = ds 2 — (i/j\H{t) 2 \i/j) dt 2 and the speed term v{s) 




FIG. 1: Generalized Fermat Principle. Natural evolution of a system from initial to final state (thick) in the projective 
Hilbert space (a), on the energy shell in classical phase space (b) and in the classical configuration space (c) together with 
fictitious paths (dashed) connecting the same states. The time functional T is an extremum along the physical path in case 
(a) and (b) but in general not in (c). The path is a geodesic on the energy shell in (b) but not in (a) and (c). 



becomes (ip\H (t)\ip) . After substituting into equation [6] 
with the new integrand and taking the inner product with 
\ip), we arrive at 

<v#) = (^> 2 + i(</#)i 2 + (m(t) 2 \*p) (9) 

We notice that both Schrodinger's equation and the non- 
linear equation [8] satisfy this equation. The extension of 
the FS metric to the new one is analogous to the tran- 
sition from the Maupertuis action J pdq to the Hamilto- 



nian action in classical mechanics J dt (pq — H(p,q,t)). 

Had we chosen the new metric to be ds' 2 = ds 2 + 

(<ip\H(t)\ipy dt 2 the new velocity term would have be- 
i 

come (ip]!!^) 2 ]^ 2 . These time dependent variational 
principles are to be contrasted with the functional 
(tp(t)\H(t) — id t \ip(t)) from which Schrodinger's equation 
can be derived. The later however lacks an appealing ge- 
ometric interpretation. 



— Fermat Principle in Phase Space We consider a 
conservative classical system composed of N particles 
described by the generalized coordinates qi,Pi, with 
velocity-independent potential V{q\, ...gjv), kinetic en- 
ergy T = ^2 i=1 p 2 and Hamiltonian H = T + V. Let the 
distance element in phase space be described by the Eu- 
clidean metric ds 2 = YliLi dq 2 -\-dp 2 . The speed along the 
physical path generated by Hamiltonian flow, v{s) equals 



Ef =1 teV 



\dq- J ^ (f?) ^ e therefore express the 

time functional T as 



T 



r x f 

J Xi 



ds 



(10) 



where Xi denotes all generalized coordinates qi^Pi treated 
on equal footing. Adding quantities of different dimen- 
sions in the metric and velocity terms should not disturb 
the reader because it simply reflects the fact that the axes 



of the phase space {qi,Pi} have different dimensions. De- 
pending on the particular problem, any of {qi,Pi} can be 
redefined in terms of the physical constants of the prob- 
lem so that they have the same dimensions. 

We investigate Fermat principle in phase space on the 
energy shell H{xi) = E (see Fig. l.b), in contrast 
to the quantum case where we didn't have to impose 
that constraint. That translates mathematically into the 
variational problem 5T = subject to the constraint 
H(xi) — E = 0. If we parametrize any arbitrary tra- 
jectory connecting and Xf by r, we can define the 
integrand L which satisfies the EL equation 



SL d SL 







(11) 



5x dr Sx 

to be L = X(r)H f (xi) + F{xi)G(xi), where H' = H-E, 
F = 



*52i ( dxi ) 



and G = x/^iX 2 . Substituting into 



equation 1 1 , we end up with 



dxi 



G 



dF 

dx. 



G 2 x\F 



V ^x 

dxj 



= 0. 



(12) 
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If Fermat principle holds in phase space, the previous 
equation has to be satisfied for every X{. By multiplying 
each equation by x\ and adding up all of them, we end 
up with 



0. 



(13) 



We see immediately that the Hamiltonian equations of 
motion q { — and pi = — |^ satisfy this equation, 
when r is taken to be the real time. As in the quan- 
tum case, the time-reversed evolution q\i = — f^- and 



Pi = will also satisfy this equation. For non-periodic 
systems, that would be of relevance to the classical ergod- 
icity problem, where the system starting from any point 
in phase space will pass arbitrarily close to any other 
point on the same energy shell via forward and back- 
ward time evolution. It is easy to check after solving the 
Hamilton's equations of motion for a given system and 



substituting back in equation [12] that we get consistent 
solutions for A(r) from all the equations. 

We point out that by the same procedure one can 
show that the phase space path length J ds is also a 
stationary quantity under the constraint of conserved 
energy for Hamiltonian flow, i.e., unlike configuration 
space, the system follows a geodesic path, at least lo- 
cally, on the enrgy shell in phase space. One interest- 
ing question to ask is whether the phase space path can 
be determined uniquely by extremizing J ds under time 
constraint rather than energy constraint. The time con- 
straint is j = T and it is an isoperimetric constraint 
rather than a holonomic constraint like the former case. 
The answer is No! One can check for any simple system 
that after substituting in the EL equation with the so- 
lution for Xi(r) generated by the Hamiltonian flow, that 
one does not get consistent solution for the Lagrange mul- 
tiplier A(t) for every X{. Therefore, counterintuitively, 
time constraint is not equivalent to energy constraint in 
phase space. 

In classical systems where kinetic and potential en- 
ergies cannot be straightforwardly defined, e.g., classi- 
cal spin systems, we hypothesize that Fermat principle 
will continue to hold in the proper space that completely 
specify the state of the system. In case of N inter- 
acting classical spins, the phase space would be a 2N- 
dimensional manifold where each spin is specified by two 
angles (^,^) and the conjugate variables {qi,Pi} would 
correspond to {cos(6i), <pi} [T2] . 

We finally comment on the validity of Fermat Prin- 
ciple in the configuration space of a conservative classi- 
cal many-body system described by the coordinates {^} 
(Fig. l.c). The metric is defined by the kinetic energy 
term which we assume to be quadratic in velocities, i.e., 
T = ^rriijqiqj. It is easy to verify that T is related to the 
speed u(s) of the Hamiltonian flow in the configuration 
space by v(s) = y/2T. The validity of Fermat Principle 



in configuration space would imply that the action J 4= 

is equivalent to the Jacobi action J ^/Tds, i.e., yields the 
same equations of motion. To test this supposition, we 
compare the Euler-Lagrange equations that will result 
from extremizing the two actions. Let ds = F (q, q) dr 
and VT = G (q, q, r). We can see by straightforward sub- 

where L is the integrand 



stitution into 

oqi 



d 5L 
dt 5q\ 



of the action, that the two actions J ^dr and J FGdr 
will satisfy the same equation of motion if F(q,q) sat- 
isfy If ~ ftWi = °' That means that / ds snould it- 
self be a stationary quantity along the Hamiltonian flow, 

or in other words the system follows a geodesic in the 
bare configuration space. Since the flow is not geodesic 
in this space except for the case of zero potential en- 
ergy [4 (or more generally when the system is confined 
to an equi-potential energy surface in the configuration 
space), we conclude that Fermat principle is generally 
not valid in this case for a non- vanishing potential , i.e., 
ST 7^ [13] . To make the previous argument more rigor- 
ous, one should take into consideration while solving the 
EL equation and comparing the two functionals the con- 
straint that the potential energy is always smaller than 
the total energy, V(qi, <?2, ••••) < E , but we suppose that 
the same conclusion will be reached. 

In conclusion, we have introduced a generalized Fer- 
mat principle as a unified concept in Hilbert and phase 
spaces of quantum and classical systems that illuminates 
new aspects of Fermat 's surmise of natural economy. 
The generalized Fermat principle provides new geometric 
variational principles that are satisfied naturally by the 
Schrodinger's and Hamilton's equations of motion. 

The concepts developed in this study provokes several 
questions which are left open. What are the implica- 
tions of the generalized Fermat principle on the nature 
of time in classical and quantum dynamics? Why is it 
natural to expect that the equations of dynamical evolu- 
tion in a complete space should yield longer times along 
all paths other than the actual path? Can this observa- 
tion be of any practical use to simulate complex many 
body systems that require unaffordable computing re- 
sources, or at least get insight into their behavior? Can 
we use the quantum Fermat principle to define a new ba- 
sis where the quantum evolution follows geodesic paths? 
What are its implications for quantum computing? Are 
there situations where the nonlinear equation of motion 
for the wave function (JsJ) that satisfies Fermat principle 
be physically relevant? What does Fermat principle say 
about the statistical evolution of an ensemble of systems 
in phase space? Can it be used to regard Liouville's equa- 
tion as a variational principle? We hope these questions 
will motivate more researchers in the future to investi- 
gate the nature of the geometry of evolution in classical 
and quantum mechanics. 
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